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Canted Ferromagnetism in Double Exchange Model with on-site 

Coulomb Repulsion 

Naoum Karchcv and Vasil Michevp 
Department of Physics, University of Sofia, 1164 Sofia, Bulgaria 

The double exchange model with on-site Coulomb repulsion is considered. Schwinger-bosons 
representation of the localized spins is used and two spin-singlet Fermion operators are introduced. 
In terms of the new Fermi fields the on-site Hund's interaction is in a diagonal form and the true 
magnons of the system are identified. The singlet fermions can be understood as electrons dressed 
by a cloud of repeatedly emitted and reabsorbed magnons. Rewritten in terms of Schwinger-bosons 
and spin-singlet Fermions the theory is [/(I) gauge invariant. We show that spontaneous breakdown 
of the gauge symmetry leads to canted ferromagnetism with on-site spins of localized and 
delocalized electrons misaligned. On-site canted phase emerges in double exchange model when 
Coulomb repulsion is large enough. The quantum phase transition between ferromagnetism and 
canted phase is studied varying the Coulomb repulsion for different values of parameters in the 
theory such as Hund's coupling and chemical potential. 
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I. INTRODUCTION 

Spin- fermion model describes materials which get their 
magnetic properties from a system of localized magnetic 
moments being coupled to conducting electrons. The 
model is known as s — d (or s~ f) model in which the elec- 
trons are separated into delocalized s electrons and local- 
ized d{f) electrons. The names of the s and d{f) elec- 
trons do not necessarily mean that the orbital electron 
states are of corresponding type. They are introduced to 
distinguish the localized from delocalized electrons. The 
model appears in the literature also as the Ferromagnetic 
Kondo Lattice model (FKLM) or the Double Exchange 
model (DEM)0,|1|1|1|1||. 

The double exchange model, as we shall call the model 
from now on, has great variety of applications to differ- 
ent topics in the magnetism. On the basis of the double 
exchange model a microscopic theory for ferromagnetic 
hexaborides, such as EuBq , is proposed ,7]. The mag- 
netism of ferromagnetic metal is found to arise from the 
half-filled 4/ shell of Eu , whose localized electrons ac- 
count for the measured moment 8, '9, 10]. The transport 
properties such as the Hall effect, magnetoresistance, and 
dc resistivity are quantitatively described within DEM 



Another materials for which the DEM is rele- 
vant are the dilute magnetic semiconductors, such as 
Gai-xMn^As, where fraction of non-magnetic elements 
Ga is replaced with magnetic ones Mn. The atoms 
of manganese supply of both carriers (holes) and mag- 
netic moment. These materials have attracted great at- 
tention after the experimental observation of ferromag- 
netic transition temperature Tc as high as llOK llj . 
One-band and two-band double exchange models have 
been solved, by means of dynamical mean-field approxi- 
mation or Monte Carlo simulations, to find the magnetic 
transition temperature as a function of coupling con- 
stants, hopping parameters, and carrier densities |l2lll3l |. 



The double exchange model is a widely used model for 
manganites [lllij. In isolation, the ions of Mn have an 
active Sd-shell with five degenerate levels. The degener- 
acy is presented due to rotational invariance within angu- 
lar momentum / = 2 subspace. The crystal environment 
results in a particular splitting of the five d-orbitals {crys- 
tal field spliting)mto two groups: the eg and t2g states. 
The electrons from the eg sector, which form a doublet, 
are removed upon hole doping. The t2g electrons, which 
form a triplet, are not affected by doping, and their pop- 
ulation remains constant. The Hund rule enforces align- 
ment of the three feg spins into a s = 3/2 state. Then, 
the t2g sector can be replaced by a localized spin at each 
manganese ion, reducing the complexity of the original 
five orbital model. The next drastic simplification is that 
only one Cg orbital is available at each site. To justify this 
one can assume 14] that a static Jahn- Teller distortion 
leads to a splitting of the degenerate Cg levels, allowing 
to keep only one active orbital. The only important in- 
teraction between the two sectors is the Hund coupling 
between localized t2g spins and mobile Cg electrons. 

The double exchange model has a rich phase diagram, 
exhibiting a variety of phases, with unusual ordering in 
the ground states. The procedures followed to obtain 
the phase diagram are different: numerical studies 15], 
dynamical mean field theory (Ifi] , and analytical calcu- 
lations [13, UMi but four phases have been systemati- 
cally observed: (i) antiferromagnetism (AF) at a density 
of mobile electrons n = 1, (ii) ferromagnetism (FM) at 
intermediate electronic densities, (iii) phase separation 
(PS) between FM and AF phases, and (iv) spin incom- 
mensurable (IC) phase at large enough Hund coupling. 
The competition between spin spiral incommensurate or- 
der or phase separation and canted ferromagnetism is also 
a topic of intensive study [VA \18L Il9| . 

The simplest but realistic Hamiltonian for the double 



exchange model has the form 






.'^ia'^j'^ 



- Jh ^ Sj • Si 



(1) 



where c^ and Ci^ are creation and destruction opera- 



tors for mobile electrons, sf = 5 X) ^ta'^a 



r'C,- 



with the 



Pauli matrices {t^jT^jT^), is the spin of the conduction 
electrons, and Si is the spin of the localized electrons. 
The sums are over all sites of a three-dimensional cubic 
lattice, and {i,j) denotes the sum over the nearest neigh- 
bors. In equation (Q the hopping amplitude and the 
Hund coupling between localized and mobile electrons 
are positive. 

The Hamiltonian (^ of the DEM is quadratic with 
respect to the fermions (c^,Cicr). Averaging in the sub- 
space of the itinerant electrons one obtains an effective 
Heisenberg like model in terms of core spins Si [20, |21j . 
In the small Jh limit Ruderman-Kittel-Kasuya-Yosida 
(RKKY) theory is recovered. The subtle point is that 
if we use a Holstein-Primakoff representation for the lo- 
calized spins Si, the creation and annihilation bose op- 
erators do not describe the true magnon of the system 
|2^. The true magnons are transversal fluctuations cor- 
responding to the total magnetization which includes 
both the spins of localized and delocalized electrons. 
Therefore the RKKY validity condition requires not only 
small Hund's coupling, but it also insists the charge car- 
rier density to be small, which in turn means that the 
magnetization of the mobile electrons is inessential. 

Since the only interaction between localized and delo- 
calized electrons is the Hund coupling, it is desirable to 
treat the on-site term in the Hamiltonian Q exactly. To 
this end, the Holstein-Primakoff transformation was gen- 
eralized to the case when the length of the spin is operator 
by itself [23. The procedure removes all spin variables 
from the Hund coupling term. The bosonic and fermionic 
sectors are constructed in - expansion up to the fourth 
order. Alternativ ely, two spin-singlet fermi fields are in- 
troduced in Ref. [2J| . In terms of the singlet Fermi fields 
the on-site term is in a diagonal form, the spin variables 
are removed, and one can treat it exactly. An analogous 
technique is used in Ref. |25l |. 

More realistic DEM would account for the on-site 
Coulomb repulsion. The Coulomb (Hubbard) interac- 
tion has a profound effect on the band structure, mag- 
netic ground state (magnetic configurations), and trans- 
port properties of the spin-fermion systems. The results 
of the electron-electron repulsion depend on parameters 
in the theory such as doping, band width, and temper- 
ature. While some effects of the Hubbard term have 
been addressed in the past by means of mean-field the- 
ory [26, 27, 28], local spin-density approximation (LSDA) 
and LSDA-fU calculations [23, and dynamical mean- 
field theory (33 its impact is not fully appreciated so 
far. 

In the present paper we study canted ferromagnetism 
in double exchange model with on-site Coulomb repul- 



sion. Usually the canted magnetism is considered as a 
two-sublattice spin configuration with neighboring lat- 
tice spins misaligned by an angle 9. First de Gennes ob- 
served that in double exchange model spin-canted state 
(Fig. ^) interpolates between ferromagnetic and antifer- 
romagnetic order [^. The canted phase (Fig. ^p) is a 
part of the phase diagram of mixed-spin (S'i)5'2) Ji — J2 
Heisenberg model on square lattice |32l |. Finally, canted 
jhase appears in the lattice models of quantum rotors 
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FIG. 1: Sketch of two sublattice spin-canted states: 

a) canted antiferromagnetism, b) canted ferromagnetism 

In the present paper we show that canted ferromag- 
netism with on-site spins of localized and delocal- 
ized electrons misaligned (Fig. ^ emerges in double 
exchange model when Coulomb repulsion is large enough. 
We study quantum phase transition between ferromag- 
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FIG. 2: Sketch of on-site spin-canted states: 
a) when Jh > 0, b) when Jh < 

netic and on-site canted orders when Coulomb repulsion 
is varied for different values of parameters in the theory 
such as Hund's coupling and chemical potential. 

The paper is organized as follows. In Sec. II we study 
double exchange model (^ supplemented with an an- 
tiferromagnetic Heisenberg interaction between nearest- 
neighbors core spins. Schwinger-bosons representation of 
the localized spins is used and two spin-singlet Fermion 



operators are introduced. In terms of the new Fermi fields 
the on-site Hund's interaction is in a diagonal form and 
the true magnons of the system are identified. The sin- 
glet fermions can be understood as electrons dressed by 
a cloud of repeatedly emitted and reabsorbed magnons. 
Integrated over the singlet fermions we obtain an effective 
Heisenberg-like theory. Positivity of the spin-stiffness, as 
a function of Hund's coupling Ju and charge carrier den- 
sity, is the condition for the stable ferromagnetism. The 
phase boundary is depicted both in the case of zero and 
non-zero antifcrromagnetic exchange. Sec. Ill is devoted 
to the on-site canted ferromagnetism in DEM with on- 
site Coulomb repulsion. The theory rewritten in terms 
of Schwinger-bosons and spin-singlet Fermions is an C/(l) 
gauge invariant theory. We show that on-site canted state 
is a state with spontaneously broken gauge symmetry. 
The quantum phase transition between ferromagnetism 
and canted phase is studied varying the parameters in 
theory. A summary in Sec. IV concludes the paper. 



II. MAGNONS IN DOUBLE EXCHANGE 
MODEL 



We consider a theory with Hamiltonian 

h = H — fJ,N = —t y^ {ctjCjcr + h.C.) — fl y^ Tlj 

(ij) i 

+ Jaf ^ S,; • Sj - Jn'^Si ■ Si (2) 

(ij) i 

where fi is the chemical potential, and Ui = cf^Cia. The 
antiferromagnetic Heisenberg term (Ja_f)0) is very im- 
portant for the manganites. In the limit when all eg 
electrons are removed, and the system is without mobile 
electrons, the t2g electrons induce an antiferromagnetic 
Heisenberg exchange between nearest-neighbors leading 
to the standard antiferromagnetism. The most promi- 
nent example is CaMnOs 14]. 

In terms of Schwinger-bosons {fi.aifla) the spin op- 
erators have the following representation: 



« 1 + 



vtaVtT = 2s. 



(3) 



The partition function can be written as a path in- 
tegral over the complex functions of the Matsub- 
ara time ficr{T) {ftdT)) and Grassmann functions 

C^.{T) (C+ (r)). 



Z(/3) = / dfi{ip+,ip,c'^,c)e 

with an action given by the expression 





-s 



(4) 



S 



dT 



T.i^tAr)^^^ir)+ct{T)c,{r)) 



h{ip+,(f,c+,c)] , 



(5) 



where /3 is the inverse temperature and the Hamiltonian 
is obtained from Eqs.Q and ^ replacing the operators 
with the functions. In terms of Schwinger-bosons the 
theory is invariant under U{1) gauge transformations 



^.W(^) 



— Jo'ji'^) 



V'jt(t); iPj^ir) =e 



= p—^°'3(^),„+ 



vUr) (6) 



with parameters which are period functions of Matsub- 
ara time aj{0) = aj{(3). The measure for the Schwinger 
bosons includes Dirac-(5 functions that enforce the con- 
straint © and the gauge-fixing condition 



D^i{^+,^) = W 



DipliT)D^,Ar) 
2m 



(7) 



\[5{vUr)^^.{r)-1s)\{5{g.J). 



We introduce two spin-singlet Fermi fields 
1 






/Ts 
1 

/2i 



¥'^(T)Cm(T), 



(8) 



[ip^l{T)c^2{T) - V3j2 (t)c,i (t)] , (9) 



which are gauge variant with charge -1 and 1 respectively 

*^^(r) = e-*"^(")*/(T), *f (r) = e^^^^'^^f (r). (10) 

The equations © and ^ can be regarded as a SU(2) 
transformation 



'^^'y =9ta'C'^'y' =^ at = 



1 



'2s 



Vtl '/'i2 



(11) 



with V&f = ^ij and *f = *2i- For that reason the 
Fermi measure is invariant under the change of variables. 
In terms of the spin-singlet Fermi fields the spin of the 
conduction electrons s, has the form 



where 






Or^^Trg+T^^g^T" 



(12) 



(13) 



It is convenient to introduce three basic vectors which 
depend on the Schwinger-bosons 



rpl r~\ /J.1 rp2 r~\ fj,2 rpZ r\ /i3 



ifj, 



tfi 



ifi 



(14) 



where Tf — -Si. Then, the spin of the electrons can be 
represented as a linear combination of three vectors Sj , 



Pj = T] + iT] and P^^ 



ls,(*f+*f-*f+*f) (15) 



The basic vectors satisiy the relations Sf = s^, Pf = 
P+2 = S, • P, = S, • P+ = 0, and P+ • P, = 2. Using 
the expression H15|l for the spin of itinerant electrons the 
total spin of the system S*°* = 8^+8^ can be written in 
the form 



1 



^+if,A 



(^qff+qi 



^B+^i 



-P,;*f+*:^ 



-P+vJ/^+vI/^ 



(16) 



The gauge invariance imposes the conditions 
(*f +*f ) = (*f +*f ) =0. As a result, the dimension- 
less magnetization per lattice site M = ((5'*°*)^) reads 



s 



i((*f+*f-*f+vl;f)) 



At zero temperature (Sf ) = s and M = s 



= ^{{^r^t 



*r*f)> 



(sn (17) 

m, where 

(18) 



is the contribution of the itinerant electrons. 

Let us average the total spin of the system (Eq. Il()l) in 
the subspace of the itinerant electrons (S*°')/ ~ M^. The 
vector Mi (M^ = M'^) identifies the local orientation 
of the total magnetization. Accounting for the gauge 
invariance, one obtains an expression for Mi in terms of 
core spins Si 



(Sf*)/ 



M, = — S, 
s 



(19) 



Now, if we use Holstein-Primakoff representation for the 
vectors M, 



Af/ = Mji + iMj2 = d2M - aja 



ijaj 



MJ 



2M ~ 



iNL 



i2 



2A/ 



a+flj 



(20) 



a^ aj 



the bose fields a,- and af are the true magnons in the 
system. In terms of the true magnons the Schwinger- 
bosons l|3l have the following representation 



An important advantage of working with A and B 
fermions is the fact that in terms of these spin-singlet 
fields the spin-fermion interaction is in a diagonal form, 
the spin variables (magnons) are removed, and one ac- 
counts for it exactly 



^s..s.=f5:[*f+*f-*f+v,f] 



(24) 



To proceed we rewrite the action (|S1) as a function of 
Schwinger-bosons and spin-singlet fermions 



S 



dr 



ip+ip,, + *f + ( ^ + ^^t^^- ] *f 



+ *■ 



d 



1 



dr 2s 



V>iaVi. 



^ 



Is 



+ /l((^+,(^,*+,*) 



(25) 



It is convenient to write the Hamiltonian h {(p'^ , (p, "^^ , '^) 
as a sum of three terms 



h = hf + IlH + hint 



(26) 



where hf is the free A and B fermions' Hamiltonian 

- ^E(*^*^-*^''*f)' (27) 



hn is the Hamiltonian of Heisenberg theory of antiferro- 
magnetism (jSJ , and hint is the Hamiltonian of magnon- 
fermion interaction 



(Pii = -\/2s- -rfajdi 



^^^-\Im 



(21) 



Replacing in Eqs. IjHJl and (|5J) for the spin-singlet 
Fermions and keeping only the first two terms in 1/Af 

l — ^^afai + ... we obtain 



expansion 



1 21/ ^i '^i 



1 



/2M 
1 

/2M 



a, ca 



— a+aiC,i + .... ,(22) 



OiCji 



^«.+ «.C,2 + 



(23) 



The equations (|22|) and H23|) show that the singlet 
fermions are electrons dressed by a virtual cloud of re- 
peatedly emitted and reabsorbed magnons. 



/lint = -*E 



^('^t'^j-- 2s) *+.*,.'+ /i.e. 



(28) 



^ (^+^+ - ^^ i^f^^f ~ *f +*f ) + h.c. 



2s 



The action (|25|l is quadratic with respect to the spin- 
singlet fermions and one can integrate them out. We can 
accomplish this by first using the representation H21I) of 
the Schwinger-bosons, then keeping only the quadratic 
terms with respect to the magnons, and finally calcu- 
lating the diagrams in the leading order of gradient ex- 
pansion. The action of the effective theory, in Gaussian 



approximation is 
ScS ^ dr 





a,- flj 



(29) 



(u) 

where M is the dimensionless magnetization per lattice 
site Eq. l(T7|) at zero temperature, and J is the effective 
exchange coupHng 





1 
1 

FM 

/ 

/ ■ 



J 



NP 



Jaf 



(30) 



71 71 71 / '^ \ 

^///(S^fe^^^'^j^ 

_Tr — -TT — -TT \A^— -L / 



-^+«f) 



— TT — TT — TT 

TT TT TT 



2i^ 



SA'PJhs 



'^''^ '^Vsin^fc 

7 Oiii fin 



(27r)3 



— TT — TT — TT 



^M=l 



In equation (p!n|l n^ — 9 {—£^) (R.=A, B) are the occupa- 
tion numbers for the A and B fermions with dispersions 



ej, = — 2i (cosfca; + cosfcj^ + cos/cx) — /i — 



e^ = —2t {cos kx + cos ky + cos kz) — i-i 



Jhs 



(31) 



Jhs 



The first term in equation H30|) comes from " tadpole" di- 
agrams with one A or B-fermion line with vertices which 
relate to the first term in the Hamiltonian of interac- 
tion H28|l . The second term is obtained calculating the 
one-loop diagrams with A and B-fermion lines, and with 
vertices which relate to the second term in hint- The 
term with time derivative in the effective action (|29|l is 
obtained summing two terms. The first one is the term 
with time derivative in the action (|25|l which in terms 

l3 

of magnons has the form J dr^a^di, while the second 



results from "tadpole" diagrams with vertices related to 
the second and third terms of the action (|25|l . 

Based on the rotational symmetry, one can supplement 
the action H29|l up to an effective Heisenberg like action, 
written in terms of the vectors Mi 



Hp 



-J^Mi-M;. 



(32) 



The ferromagnetic phase is stable if the effective ex- 
change coupling constant is positive J)0. The dimension- 
less constant J/W, where W — 12i is the band-width, de- 
pends on Jaf/W, Jhs/W and iJ./2t. The (^^,n) phase 
diagram, where n is carrier density, is depicted in Fig. 3 
for Jaf = and ^ = 0.1. 

The phase diagram Fig. 3a [Jaf = 0) is in a good 
agreement with phase diagrams obtained numerically [l5j 



FIG. 3: Phase diagrams when a) Jaf = and b) -^ — 0.1 



and by means of alternative analytical calculations [T^ ■ 
Phase diagram Fig. 3b shows that direct antiferromag- 
netic exchange suppresses the ferromagnetism at small 
values of carrier concentrations, which is well known ex- 
perimental fact for manganites \i^ . 



III. CANTED FERROMAGNETISM 

After considering pure double exchange model, let us 
address the double exchange model supplemented with 
on-site Coulomb repulsion (Hubbard term). 

h = H - fiN = -t^{c+^cj„ + h.c.) 

- JH^Si-s, + U^ni^nii-fj.^n„ (33) 



where riicr = c^^Ci^- Our purpose is to show that canted 
ferromagnetism, with on-site spins of localized and delo- 
calized electrons misaligned, emerges in double exchange 
model when Coulomb repulsion is large enough. 

Let us average the spin of the electrons (Eq. I15|l in 
the subspace of the itinerant electrons. We obtain, as a 
consequence of gauge invariance, that spin of electrons 
are parallel to the localized spins (s^)/ = — S.^. Equa- 
tion IjlSI) shows that the on-site spins are misaligned if 
(vj/^+vI/S^ and (vl/f +vl/^) are not equal to zero, which 
in turn means that gauge symmetry is spontaneously 
broken. To explore dynamical breakdown of the gauge 
symmetry, or which is the same, the on-site canted ferro- 
magnetism we rewrite the Hamiltonian (I33|l in terms of 
Schwinger-bosons and spin-singlet Fermions. In particu- 
lar, one obtains for the Hubbard term 



E 



n,Tnu = -E*^*f*^*^ 



(34) 



We decouple this term by means of the Hubbard- 
Stratanovich transformation, introducing complex field 



Ai (A^), the order parameter of the gauge symmetry 
breaking. 



e ' 

r '^ 

= dfi{A+A)exp ~ / '^'^X! 



(35) 

A+(r)A.(r) 
U 



+ *f+(T)vI/f (r)A,(r) + A+(r)*f+(r)vI/f (r) 



Now, the partition function |0J can be represented as a 
path integral over the spin-singlet fermions, Schwinger- 
bosons, and complex order parameter. The integral over 
the fermions is Gaussian, and one can integrate them out. 
The resulting expression for the partition function is an 
integral over the Schwinger-bosons, and complex order 
parameter. 

Z(/3) = /d^(^+,^,A+,A)e-^(^^'^'^^'^). (36) 



We perform the integral over the collective variables A^ 
and Ai using the steepest descend method. To this end, 
we set the spin fluctuation af and Ui (see equations 
(|20I21|| ) equal to zero and assume that mean- field value 
of the order parameter Aj(r) is an independent of r and 
lattice sites i real constant A. Then, the free energy of 
the system, in mean-field approximation, is 



T 



U 



-T 



f 



(37) 



where J^f is the free energy of a Fermi system with Haniil- 
tonian 



h = J2 £^*r*nef*f+vi,f (38) 

k - 

To write the Hamiltonian in diagonal form one introduces 
new Fermi fields ^^ and tp'l 



^^^urk + vi'l 

where the coefficients are 



*f 



^vrk + ui>l (39) 



1 



Jh s 



\2\ ^{JHsf+iA^J' 



Then, 



(signA)-\/l — V? 



hf^T.H'('k''^k+^li'l^i'l] 



(40) 



(41) 



Here ef. = e,. , st — et , where 



fc ' 



e^^ek±;;\/{JHsr+4A^, 



(42) 



and Ek = —2t [cos(fca;) -I- cos{ky) + coa{kz)] — fj.. Now, we 
can obtain the mean-field expression for the free energy. 
At zero temperature it is 






TT 7T 7T 



(27r)3 



[eie{-el)+el0{-el)] (43) 



— 7r — 7r — 7r 



It is convenient to introduce the angle between the on- 
site spin of carrier and localized spin 



COS0 ~ 



^i ' S; 



(44) 



Using the equation (|15|l . for the spin of itinerant elec- 
trons, and the properties of the basic vectors one obtains 



COS0 — 



4(vI/f+M.f)(vI/f+*^^ 



(*f +^^ - *f +^f ) 



(45) 



We calculate the matrix elements in the formula H45|) 
in mean-field approximation applying the transformation 
(EnEDJ). The result is 



(^i4+^A_ ^S + ^S\ 



l^f+^f) = i^f+^f) = 



Jh s {n" - n^) 
[Jh sf + 4A2 
A fn" - n'>) 



(46) 



(Jhs) 



-4A2 



where n" and n^ are occupation numbers for " a" and " b" 
fermions introduced by the transformation 139(1 . After 
some algebra we arrive at the mean-field expression for 
the angle 



cos 9 



Jhs 



[Jh sy + 4A2 



(47) 



Next, we replace A in equations H42|l and (|43|l by cos0 
from fEa l47|) and rewrite the mean-field free energy as a 
function of cos 8. The dimensionless energy F = 6J-'/W 
is depicted in (Figs .4, 5, 6) for different values of W/U 
and fixed Jh s/W and n/W. As graphs show, increasing 
the Coulomb repulsion constant the system passes trough 
a first order quantum phase transition. Red lines corre- 
spond to the critical values Uc of the Coulomb repulsion. 
The values Uc and Oc depend on the parameters of the 
theory such as Hund's coupling constant, chemical po- 
tential and band width . The character of the transition 
also depends on the parameters in the theory. We see 
(Fig. 6) that when Jh s/W = 0.50 and fi/W = -0.33 
for small values of Coulomb repulsion, W/U=1.50, 1.10, 
the minimum of the free energy is at cos = 1. Near the 



quantum phase transition W/U = 0.77 the ground state 
is highly degenerated, while below this critical value, for 
large enough U, the on-site canted ferromagnetic state, 
with cosQ < 1, is the ground state of the system. To 
figure that quantum phase transition out one has to go 
beyond the mean-field theory or to use alternative meth- 
ods of calculations to complement the above one. 




Cose 

FIG. 4: Dimensionless mean-field free energy F - 
a function of cos0 for Jh s/W = 2.32, fi/W = 
W/U^O.23; 0.25; 0.28; 0.30; 0.32; 0.35; 0.40 



6T/W as 
-1.22, and 




FIG. 5: Dimensionless mean-field free energy F — 6J-/W 
as a function of cos 6* for Jhs/W = 0.95, fJ^/W = —5, and 
W/U^0.30; 0.33; 0.36; 0.39; 0.42; 0.45; 0.48 




Cose 

FIG. 6: Dimensionless mean-field free energy F = 6J-/W as 
a function of cos 6 for Jh s/W = 0.50, ^J./W — —0.33, and 
VF/f/=0.32; 0.39; 0.48; 0.60; 0.77; 1.10; 1.50 



IV. CONCLUSIONS 

We have argued that the on-site Coulomb repulsion 
strongly affected the magnetic properties of spin-fermion 
systems. In particular, when Coulomb repulsion is strong 
enough, the on-site localized and carriers' spins become 
misaligned (on-site canted ferromagnetic state). As fol- 
lows from (EqEI), cos 6 > when J^r > (see Fig. 
2a), and cos 9 < when Jat < (Fig. 2b), To ob- 
tain this result a double exchange model with Hubbard 
term was considered. We represented the localized spins 
by means of Schwinger-bosons and introduced two spin- 
singlet Fermion operators. In terms of the new Fermi 
fields the on-site Hund's interaction is in a diagonal form 
and the true magnons of the system can be recognized. 
Written in terms of Schwinger-bosons and spin-singlet 
fermions the theory is U{1) gauge invariant. We have 
shown that on-site canted ferromagnetic state is a state 
with spontaneously broken gauge symmetry. This is be- 
cause the order parameter is gauge varying collective field 
with charge -2 (see equations (|10() and (|35|l ') 



A'(r) 



'MA,(r) 



(48) 



and non-zero expectation value (Aj(t)) ^ means spon- 
taneous breakdown of the gauge symmetry. 

To study the Goldstone modes in the on-site canted 
ferromagnetic phase it is convenient to represent the 
Schwinger-bosons and the order parameter in the form 



s''^^/2s 



Vjl 



A, = |A,|e''^^-. 



M 



a J a-i, 



(49) 
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Then the new fields, af , Uj, (f>j and Xj transform under 
the gauge transformation in the following way 



xj, 



JCj, 



X, 



Xj - 2aj 



(50) 



The theory is an C/(l) gauge theory and we have to im- 
pose one gauge fixing condition. Hence, there are two 
Goldstone modes in the theory. For example one can 
use (j)i = as a gauge fixing condition (see equation 
(|21|l ). then the Goldstone modes are the magnons at (a^) 
and Xi phase. The physical origin of the extra mode is 
the totaly broken rotation symmetry, while mathemat- 



ical reason is the spontaneous breakdown of the gauge 
symmetry. Alternatively, one can choose Xi — for the 
gauge fixing. Both these conditions are not convenient. 
In the quadratic parts of the corresponding effective the- 
ories there are terms which mix magnons and phases. 
These terms are obstacle to recognize the spectrum in the 
theory. One hopes that there is a gauge fixing condition 
which involves all gauge varying fields and the quadratic 
terms of the effective theory takes diagonal form. This 
issue will be addressed elsewhere. 
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